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Abstract 

Although self-incompatibility (SI) in plants has been studied extensively, far less is known about interspecific repro- 
ductive barriers. One interspecific barrier, known as unilateral incongruity or incompatibility (Ul), occurs when species 
display unidirectional compatibility in interspecific crosses. In the wild tomato species Solanum pennellii, both SI and 
self-compatible (SC) populations express Ul when crossed with domesticated tomato, offering a useful model system 
to dissect the molecular mechanisms involved in reproductive barriers. In this study, the timing of reproductive barrier 
establishment during pistil development was determined in SI and SC accessions of S. pennellii using a semi-Zn vivo 
system to track pollen-tube growth in developing styles. Both SI and Ul barriers were absent in styles 5 days prior to 
flower opening, but were established by 2 days before flower opening, with partial barriers detected during a transi- 
tion period 3-4 days before flower opening. The developmental expression dynamics of known SI factors, S-RNases 
and HT proteins, was also examined. The accumulation of HT-A protein coincided temporally and spatially with Ul 
barriers in developing pistils. Proteomic analysis of stigma/styles from key developmental stages showed a switch in 
protein profiles from cell-division-associated proteins in immature stigma/styles to a set of proteins in mature stigma/ 
styles that included S-RNases, HT-A protein and proteins associated with cell-wall loosening and defense responses, 
which could be involved in pollen-pistil interactions. Other prominent proteins in mature stigma/styles were those 
involved in lipid metabolism, consistent with the accumulation of lipid-rich material during pistil maturation. 

Key words: pistil development, reproductive barriers, self-incompatibility, Solanum pennellii, stigma/stylar proteins, unilateral 
incongruity/incompatibility. 



Introduction 

Flowering plants rely predominantly on generalist insect transfer of pollen from one individual to another (Faegri and 
pollinators or other non-specific vectors, such as wind, for the van der Pijl, 1979; Waser and Ollerton, 2006). Consequently, 

Abbreviations: ABF, Aniline Blue fluoroehrome; ANOVA, analysis of variance; CI, confidence index; CRP, cysteine-rich protein; iTRAQ, isobaric tag for relative and 
absolute quantitation; LTP, lipid-transfer protein; PVDF, polyvinylidene diflouride; QTL, quantitative trait locus; SC, self-compatible/self-compatibility; SDS/PAGE, 
sodium dodecyl sulphate/polyacrylamide-gel electrophoresis; SI, self-incompatible/self-incompatibility; SP, signal peptide; Ul, unilateral incongruity/incompatibility. 
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they have evolved both prezygotic and postzygotic reproduc- 
tive barriers to prevent undesirable genomic unions. Prezygotic 
barriers may contribute more than other mechanisms to out- 
crossing or reproductive isolation in plants (Rieseberg and 
Willis, 2007; Lowry et al, 2008). In the Lycopersicon section 
of Solarium, prepollination prezygotic barriers can limit the 
transfer of pollen from one species to stigmas of other species 
by mechanisms such as pollinator preference and floral mor- 
phology traits (Rick et al, 1978), while postpollination prezy- 
gotic barriers act after pollination but before fertilization. 

Self-incompatibility (SI) is a postpollination prezygotic 
intraspecific barrier, wherein self pollen is actively rejected by 
stigmas or styles to prevent inbreeding depression and facili- 
tate outcrossing. SI systems are well studied and these studies 
have provided many insights into the mechanisms that pre- 
vent self-pollination (e.g., Stephenson et al., 2000; McCubbin 
and Kao, 2000; Kachroo et al., 2002; Franklin-Tong, 2008; 
McClure, 2009; Castric et al., 2010). In S-RNase-dependent 
gametophytic SI prevalent in the Solanaceae, cytotoxic 
S-RNases and the asparagine-rich HT proteins in styles play 
key roles in pollen rejection (McClure et al., 1999; O'Brien 
et al., 2002; McClure and Franklin-Tong, 2006). Unilateral 
incongruity or incompatibility (UI) is an interspecific barrier 
that results in rejection of genetically distant pollen genotypes 
and thereby prevents genomic unions that may undermine 
fitness. As the name indicates, pollen rejection in UI is one- 
sided. In most cases, pollen tubes of SI species reach the ovary 
of SC species, while in reciprocal crosses, pollen tubes of SC 
species stop growing prior to reaching the ovary of SI species 
(Lewis and Crowe, 1958). This general trend, called the SI x 
SC rule, as well as results from genetic and molecular/trans- 
genic analysis, led to the suggestion that UI may simply be an 
elaboration of SI (Pandey, 1981; Chetelat and DeVerna, 1991; 
Murfett et al., 1996). However, the role of S-RNases in UI is 
complex, in that there appear to be both S-RNase-dependent 
and S-RNase-independent mechanisms for interspecific pol- 
len rejection (Murfett et al., 1996, Covey et al, 2010). In 
Nicotiana, where only a single HT protein is known, this pro- 
tein is implicated in UI (Hancock et al., 2005). Solarium spe- 
cies have at least two highly similar paralogous genes, HT-A 
and HT-B (O'Brien et al, 2002; Covey et al, 2010). In the 
tomato clade, these genes are tandemly duplicated at a locus 
on chromosome 12 that is linked to a UI quantitative trait 
locus (QTL; Bernacchi and Tanksley, 1997; Covey et al. , 2010). 
It should be noted, however, there are many exceptions to the 
SI x SC rule, where SC species, or SC populations within a 
predominantly SI species, exhibit interspecific pollen rejec- 
tion (Camadro and Peloquin, 1981; Liedl et al, 1996; Covey 
et al, 2010) or SI taxa accept pollen from SC relatives (de 
Nettancourt, 2001). Previous work showed that pollen tube 
growth differed in SI and UI in Solarium pennellii, and that 
a SC accession of S. pennellii could still reject interspecific 
pollen tubes (Liedl et al, 1996). As early as the 1960s, Grun 
and Aubertin (1966) surmised from their analysis of UI in 
Solanum species that genes other than SI alleles may contrib- 
ute to UI. Thus, the relationship between SI and UI remains 
unclear based on genetic data (Chetelat and DeVerna, 1991; 
Bernacchi and Tanksley, 1997; Li et al, 2010 and references 



therein). Moreover, there are no compelling data indicating 
that all instances of UI are mechanistically related, and evi- 
dence suggests that multiple mechanisms can contribute to UI 
even between a single pair of species (Murfett et al, 1996). 

We have initiated cellular and molecular characterization 
of interspecific prezygotic mating barriers in the tomato 
clade (Bedinger et al, 2011; Covey et al, 2010). Section 
Lycopersicon of the genus Solanum is particularly well suited 
for this analysis. It is a relatively small monophyletic clade 
consisting of 13 closely related diploid species (domesticated 
tomato Solanum lycopersicum and 12 wild species) with varia- 
ble mating systems both within and between species (Spooner 
et al, 2005; Bedinger et al, 2011). 

Here, we examined the developmental dynamics of UI 
and SI in pistils in three accessions of Solanum pennellii. 
In addition, developmental profiles of proteins implicated 
in SI (S-RNases and the asparagine-rich HT proteins) were 
assessed using immunological analyses. Furthermore, prot- 
eomic profiling of stigma/styles from pre-barrier and post- 
barrier developmental stages [i.e., 5 days (-5) and 1 day (-1) 
before flower opening, respectively], validated our immuno- 
logical analyses of SI factors in developing pistils, and identi- 
fied a substantial population of proteins present in immature 
and mature stigma/styles. 

Materials and methods 

Plant material and growth conditions 

One SC (LA0716) and two SI (LA2560 and LA1340) accessions of 
S. pennellii were selected for study. Domesticated tomato S. lycoper- 
sicum (cv. VF36 or M82) was used as the pollen donor in UI crosses. 
Plants were grown in Fort Collins, CO, USA, in either in the field 
(summer 2007) or in the greenhouse (2007-2010) in 25cm pots filled 
with Pro-Mix BX, watered regularly and supplemented weekly with 
Peter's 10-20-10 fertilizer. Natural day-length conditions in the green- 
house were supplemented to a 16h photoperiod. Average tempera- 
tures ranged from 27 to 32 °C during the day and were °21 °C at night. 
In the field, pollinators were excluded by bagging emasculated flowers. 

For proteomic analysis, plants were grown both in Colorado as 
described above and in greenhouses in Ithaca, NY, USA, in order to 
generate sufficient material. No phenotypic differences were observed 
between the plants grown in the two locations. Stigma/styles collected 
in both locations were pooled to minimize the biological variation 
related to the different locations. Unpollinated stigma/styles from 
S. pennellii SI LA2560 at two different developmental stages (1 and 
5 days before the onset of an thesis; termed stages -1 and -5, respec- 
tively) were collected and frozen at -80 °C until processing. 

Chemicals 

All chemicals, unless otherwise stated, were obtained from Sigma- 
Aldrich (St Louis, MO, USA). 

Semi-\n vivo pollinations and style measurements 

A semi-/« vivo system was developed to monitor pollen-tube growth 
in excised pistils laid on a solid medium [4% polyethylene glycol 
8000, 1.6mM boric acid, 4% sucrose, 3mM Ca(N0 3 ) 2 .4H 7 0, 0.8mM 
MgS0 4 .7H 2 0, ImM KNO s , 20mM MES, pH 6.0 and 2% agarose] 
in a petri dish (Fig. IB). Inflorescences with developing flower buds 
at stages -1 to -5 (days before the onset of anthesis) were used for 
pollinations («=48 inflorescences from a minimum of 3-5 plants per 
accession). Petals and sepals were removed, flowers were emasculated 
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and stigmas were pollinated with appropriate genotypes. Pistils 
were placed on the medium immediately (within seconds) after their 
stigmas were dipped in pollen. Additionally, a paint brush dipped 
in pollen was used to ensure that stigmas of all pistils in the plate 
were loaded with pollen. Plates were incubated in a humid chamber 
at 22 ±2 °C for 24 h. Sterilization of tissue was deemed unnecessary 
since no significant microbial growth was observed within this time 
period. Pollinated pistils were collected after 24 h and fixed, cleared 
and stained using Aniline Blue fluorochrome {ABF: 4',4-[car- 
bonyl bis(benzene-4, l-diyl)-bis(imino)]-bisbenzensulphonic acid; 
Biosupplies, Bundoora, VIC, Australia} in 0.1 M K 2 P0 4 as a stain for 
callose in pollen-tube walls (Kearns and Inouye, 1993) as described in 
Covey et al. (2010). Slide-mounted pistils were examined using a Leica 
DM 5500 B microscope (Leica Microsystems, Buffalo Grove, IL, 
USA) running IPLab version 4 software (BD Biosciences, San Jose, 
CA, USA) coupled with a C4742-95 camera (Hamamatsu Photonics, 
Hamamatsu-City, Japan). Images were captured at 100x magnifica- 
tion using a DAPI filter cube (425 nm long pass; Leica Microsystems). 
For each pistil (SIn=61; UIn=136), style length was measured from 
the base of the style to the top of the stigma surface. In each style, 
the 10 longest pollen tubes were measured from the end of the pollen 
tube to the top of the stigma and a mean was calculated for use in our 
analyses (e.g. mean pollen-tube length). When pollen tubes reached 
the ovary, their length was considered equal to the style length. In a 
few cases fewer than 10 pollen tubes could be measured, in which case 
we calculated the mean of the measureable tubes. All style and pollen- 
tube lengths were measured using the programme MicroMeasure ver- 
sion 3.3 (http://rydberg.biology.colostate.edu/MicroMeasure). 

We compared mean style length and mean pollen-tube length at 
five bud stages in SI and UI pollinations with four separate analyses 
of variance (ANOVAs). For each model we included bud stage (-5 
to -1), accession ID (SI, LA2650 and LAI 340; UI, LA2650, LAI 340 
and LA0716) and the interaction term as main effects using the aov 
(ANOVA) function in the R statistical package (R Development Core 
Team, 2010). In all four models, there was no significant interaction 
between bud stage and accession ID, thus the interaction term was 
dropped. Post hoc comparisons of means were made using Tukey's hon- 
estly significant difference (HSD; R Development Core Team, 2010). 



On-plant pollinations 

Pollinations were carried out with flowers still on the plant (UI n=18 
pollinations; SI «=40 pollinations). Unlike the procedure used for semi- 
in vivo studies, only anthers were removed, leaving sepals and petals 
intact on the flower in order to limit desiccation and mechanical damage. 
Trilinolein (50mg), a C 18 lipid known to promote pollen germination 
(Wolters-Arts et al., 1998), was emulsified in 0.5ml of 14% sucrose 
and 0.01% boric acid (±3mM calcium nitrate) and used to compen- 
sate for the lack of exudate on stigmas of young pistils. Stigmas were 
painted with this emulsion using a camel-hair brush before applying 
appropriate pollen samples. Pistils were collected after 24 h and fixed 
immediately for subsequent processing and staining as described in 
Covey et al. (2010). 



Microscopy of pistil tissue 

To examine pistil structure at different developmental stages, devel- 
opmental series of SI LA2560 pistils were fixed in glutaraldehyde, 
post-fixed and stained with osmium tetroxide, dehydrated with etha- 
nol and embedded in Eponate 12 (Ted Pella, Reading, CA, USA). 
Sections of 3 um thickness were cut using an ultramicrotome with a 
diamond knife. Sections presented in the manuscript were from the 
midsagittal plane of the stigmas. No additional stains were used and 
images were collected using a Leica DM5500 B microscope. 

Immunoblot analysis 

Stylar proteins were extracted and used in immunoblot detection 
of S-RNase Pen-1, HT-A and -HT-B, as described in Covey et al. 



(2010). In brief, stigma/styles were homogenized in 2 x sodium dode- 
cyl sulphate/polyacrylamide-gel electrophoresis (SDS/PAGE) sample 
buffer, boiled for 5min and centrifuged (lOmin, 14000g). Proteins 
were separated in a 15% Tris/Tricine gel, blotted onto polyvinylidene 
diflouride (PVDF) membranes and immunostained. Rabbit Anti- 
Spen-i-RNase antibody was raised against the synthetic peptide 
Ac-CTPSKTNTVIDVEKAIKAV-amide, found in S pen . r RNase 
from SI S. pennellii LA2560 (Covey et al, 2010). The affinity-puri- 
fied anti-peptide antibody was used at 1:10000 dilution. HT proteins 
were extracted and detected as described in Covey et al. (2010). 

For the immunoblot in Fig. 6, proteins were separated in a 12.5% 
Tris/Tricine gel, blotted onto PVDF membranes and immunostained 
with rabbit anti-S-RNase C2 domain antibody. This antibody was 
raised against the synthetic peptide Ac-NFTIHGLWPD-amide 
that corresponds to the conserved C2 domain of S-RNases (Haring 
et al, 1990), and the antibody was affinity purified and used at 
1:5000 dilution. 



Proteomic analysis 

Comparative isobaric tag for relative and absolute quantitation 
(iTRAQ) proteomic analyses of unpollinated stigma/styles from 
S. pennellii LA2560 at the -5 (UI and SI barriers absent) and -1 
(UI and SI barriers present) developmental stages was conducted 
to assess molecular changes during pistil maturation. iTRAQ is a 
gel-free comparative proteomic platform that allows the simultane- 
ous identification and relative quantitation of proteins from multi- 
ple samples based on the use of unique isobaric tags that identify 
the sample from which a protein originated (Wise et al, 2007). 
Proteins were extracted from three equal aliquots of pooled stigma/ 
styles as previously described (Covey et al, 2010) and then reduced, 
alkylated, proteolytically digested with trypsin and labelled with 
iTRAQ-4plex reagents, according to the manufacturer's instruc- 
tions (Applied Biosystems, Foster City, CA, USA; www.applied- 
biosystems.com/). After labelling, samples were then fractionated 
by strong cation exchange (SCX) chromatography using an Agilent 
1 100 HPLC with a UV detector (Agilent Technologies, Santa Clara, 
CA, USA). The tryptic peptides were dissolved in buffer A (10 mM 
potassium phosphate, pH 3.0, 25% acetonitrile), prior to SCX sepa- 
ration using a PolyLC PolySulfoethyl A column (2.1 mm x 150 mm; 
PolyLC, Columbia, MD, USA) eluted with lOmM potassium 
phosphate (pH 3.0), 25% acetonitrile with a 0-500 mM KC1 gradi- 
ent. During this elution 40 fractions were collected at a flow rate of 
200 ul/min. These were then pooled to create 20 uniformly complex 
second-dimension fractions based on the UV trace at 214nm. After 
solid-phase extraction (Waters SepPak C 18 cartridge; 1 ml 75% ace- 
tonitrile eluent) each fraction was then further fractionated by nano- 
scale reverse-phase chromatography on an LC Packings UltiMate 
integrated capillary HPLC with a Dionex (Dionex, Sunnyvale, CA, 
USA) Switchos valve switching unit and a PepMap C 18 RP nanocol- 
umn coupled to an LTQ Orbitrap Velos hybrid FT mass spectrome- 
ter system (Thermo Scientific, Rockford, IL, USA). This system was 
used for peptide identification with quality parameters including a 
95% confidence index (CI) cutoff, a false discovery rate of 0.05, a P 
value of 0.01 and at least two peptides per protein. Protein extraction 
and iTRAQ analysis was performed three times with equal aliquots 
of pooled stigma/styles at each developmental stage. For each of the 
three iTRAQ experiments about 100 stigma/styles at the -1 stage 
and about 200 stigma/styles at the -5 stage were used, and extracts 
were divided into two replicates for each of the three experiments. 
Those proteins showing more than a 2-fold difference in abundance 
in at least two iTRAQ experiments were considered to be differen- 
tially expressed. The resulting mass spectra data from the iTRAQ 
runs, corresponding to the different peptides, were analysed using 
Mascot 2.2 (www.matrixscience.com/) with a protein database built 
using transcriptomic data derived from 454-based cDNA sequencing 
(RNA-Seq) from several tomato tissues (stigma/styles, pollen, leaves 
and fruit; http://ted.bti.cornell.edu/cgi-bin/TFGD/digital/experi- 
ment.cgi?ID=D002) together with a customized in-house database 
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comprising S-RNase sequences annotated in GenBank. Proteomic 
data can be accessed at www.irbtomato.org. Protein identity was 
determined using both GenBank (http://blast.ncbi.nlm.nih.gov/ 
Blast. cgi) and the tomato genome sequence (The Tomato Genome 
Consortium, 2012). 

N-terminal signal peptide predictions were made using SignalP 4.0 
software (www.cbs.dtu.dk/services/SignalP/). For these analyses, pre- 
dicted protein sequences were from the best matched, most complete 
codon sequences in GenBank (Supplementary Tables SI and S2). 

Results 

Mating systems and developmental staging of 
S. pennellii 

S. pennellii, a member of section Lycopersicon - the tomato 
clade - has both SI populations, represented here by acces- 
sions LA2560 and LA1340, and SC populations, represented 
by accession LA0716. S. pennellii, therefore, provides a useful 
platform to determine the interrelations and molecular mech- 
anisms regulating both SI and UI reproductive barriers. All 
accessions of S. pennellii used in this study are interfertile and 
all display UI, in that they reject pollen from SC cultivated 
tomato, S. lycopersicum, while the reciprocal pollination is 
compatible (Covey et al, 2010). 

Our observations of both SI LA2560 and SI LA1340 show 
a consistent pre-anthesis floral developmental pattern with 
buds of successive nodes on a single inflorescence being one 
day apart in development (Fig. 1A). The node with a freshly 
opened flower is considered as stage 0 and the subtending 
nodes are successively staged as -1, -2, -3, etc. In SC acces- 
sion LA0716, petals begin to separate earlier in bud develop- 
ment, so in this accession stage 0 was assigned to flowers when 
petals (and anthers) attain full bright yellow color and were 
completely open. Gradziel and Robinson (1989) also used 
nodal position as an indicator of bud development in Solanum 
peruvianum but did not correlate this with the chronological 
age of the buds. Style length of fresh pistils was measured 
at each stage (-5 to -1) for each accession (Supplementary 
Fig. SI). As expected, all three accessions showed an increase 
in style length with bud stage (age), and the growth rate as 
determined by the slope of the relationship between develop- 
mental stage and style length ranged between 1.1 and 1.3 mm/ 
day for the three different accessions throughout the 5-day 
developmental period under study. 

Temporal and spatial onset of reproductive barriers in 
late pistil development 




Fig. 1. Developmental staging and semi-/n vivo pollination system. 
(A) Developmental staging of flower buds in an inflorescence of 
S. pennellii SI LA2560. Bud break or stage 0 indicates a flower on 
the first day of its opening. Successive buds are chronologically 
separated from each other by one day. (B) Pollinated pistils 
incubated on the medium used in the semi-/n vivo assay for pollen- 
tube growth. 



Using a semi-/« vivo approach to facilitate pollen adherence 
and germination at early bud stages (Fig. IB), we tracked the 

appearance of SI and UI reproductive barriers during pistil of -3 or younger stages to improve pollen retention and ger- 

development by pollinating pistils at different developmental mination. However, even with trilinolein application many 

stages and observing pollen tubes using ABF staining in fixed on-plant pollinations of young pistils demonstrated highly 

pistils (Fig. 2A and 2B). A limited number of on-plant polli- variable pollen adhesion, germination and pollen-tube 

nations were performed to confirm the findings with the semi- growth, highlighting the value of a reliable semi-/« vivo assay 

in vivo system (Fig. 2C and 2D). Trilinolein, a major lipid for developmental analyses. 

component of stigma exudates in Solanaceae (Matsuzaki In both SI and UI pollinations, successfully germinated 

et al., 1983; Wolters-Arts et al., 1998), was applied to stigmas pollen tubes were found to reach the ovary at early bud stages 
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Semi-/V7 vivo SI 



1 mm 




Fig. 2. Pollen-tube growth in Ul and SI pollinations. (A) S. pennellii accession SI LA2560 pollinated with S. lycopersicum cultivar M82 
pollen (Ul) in the semi-/n vivo system. (B) S. pennellii accession SI LA1340 self pollinated (SI) in the semi-/n vivo system. (C) S. pennellii 
accession SC U\0716 pollinated by S. lycopersicum cultivar M82 (Ul) on-plant. (D) S. pennellii accession SI LA1340 self-pollinated (SI) 
on-plant. All pistils were stained with ABF. Arrowheads indicate the point at which most pollen tubes have stopped growing, and arrows 
indicate the position of the longest pollen tube. VB, vascular bundles. 



(-5, -4, -3), and failed to reach the ovary at later bud stages 
(-1, -2) in the semi- in vivo system, suggesting that reproduc- 
tive barriers are not fully in place 3-5 days before flowering 
(Fig. 3A and 3B). A similar pattern of pollen-tube growth 
was found in on-plant Ul and SI pollinations, with the excep- 
tion of bud stage -3, which is likely a transition period for the 
onset of reproductive barriers (Fig. 3C and 3D). In compat- 
ible pollinations (SC LA0716 self-pollination and SI LA2560 
or LAI 340 sib pollination), pollen tubes reached the ovary at 
all developmental stages (data not shown). 

In a more detailed analysis of pollen-tube length relative 
to style length, no differences were found among accessions 
in either style length or mean pollen-tube length by bud stage 
(Table 1). In Ul crosses, mean pollen-tube length (calculated 
as the mean of the 10 longest pollen tubes/pistil) was signifi- 
cantly shorter at stages - 1 and -2 compared to earlier stages 



(Fig. 4A, Table 1), suggesting complete onset of Ul barri- 
ers by bud stage -2. Style length was significantly longer in 
later stages compared to earlier stages, but this is not likely 
to be an important factor contributing to reproductive isola- 
tion in these Ul crosses, because the net length reached by 
pollen tubes in Ul rejection is less than the length of a -5 
stage style. 

In SI pollinations, mean pollen-tube length was more vari- 
able than in Ul crosses, and the timing of the onset of SI 
reproductive barriers was not reflected in an abrupt change 
in mean pollen-tube length (Fig. 4B, Table 1). As in the Ul 
crosses, style length was significantly longer in later stages (-1 
and -2) relative to early stages (Fig. 4B, Table 1). In SI, style 
length in conjunction with the known developmental regula- 
tion of SI factors (McClure et al, 1999 and see below) prob- 
ably contributes to the onset of SI reproductive barriers. 
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Fig. 3. Pollen tubes at ovaries in Ul and SI pollinations of developmentally staged pistils. The percentage of pollinations that resulted 
in pollen tubes reaching the ovary (black) or with no pollen tubes reaching the ovary (white) at five bud stages in semi-//? vivo Ul (A) and 
SI (B) pollinations and in on-plant Ul (C) and SI (D) pollinations. The number of images analysed for each type of pollination is shown in 
parentheses. 



The location in the style of reproductive isolating barriers 
differed between Ul and SI, as has previously been reported 
(Hardon, 1967; Gradziel and Robinson, 1989; Kuboyama 
et al, 1994; Liedl et al, 1996; Covey et al, 2010) (Fig. 4). At 
later stages (-1 and -2) in SI pollinations, pollen-tube growth 
ceased in the mid-region of the style, whereas in Ul crosses 
pollen stopped in the upper third of the style. 

SI factor accumulation during pistil development 

S-RNases and HT proteins are essential SI factors produced in 
the style. We prepared antibodies specific for HT-A, HT-B and 
an S-RNase previously identified in LA2560, Pen-1 (Covey 
et al, 2010), to enable comparison of the expression profile 
of these proteins with the development of pollination barriers 
(Fig. 5). Our results show that considerable S-RNase protein 
was present in stigma/styles of SI S. pennellii LA2560 at stage 
-4 (Fig. 5A), suggesting that S-RNase alone is not sufficient 
for full rejection of either self or interspecific pollen, since 



barriers are not fully established until stage -2 (Figs 2-4). In 
fact, S-RNases are not detected in the SC accession LA0716 
(Fig. 6), which exhibits Ul pollen-tube rejection as rapidly as 
the SI accessions (Liedl et al, 1996; Covey et al, 2010). In both 
SI and SC S. pennellii accessions, HT-A protein was detected 
from stage -3 onwards, while HT-B was not detectable until 
stage -1 (Fig. 5). Although proteomic analysis (presented 
below) confirms the low abundance of HT-B relative to HT-A, 
the apparent temporal difference in HT-A and HT-B expression 
seen here could also be attributed to a relatively lower sensitiv- 
ity of the HT-B antibodies compared to the HT-A antibodies 
used in this study (i.e. trace levels of HT-B present at earlier 
stages may not have been detected). In summary, both SI and 
Ul pollen rejection commences by stage -2, when significant 
levels of both HT-A and S-RNase proteins were detected. 

Fig. 7 compares HT-A protein accumulation in upper and 
lower portions of the stigma/style in S. pennellii accessions SI 
LA2560 and SC LA0716. The higher abundance of HT pro- 
tein in the upper half of the stigma/style corresponds to the 
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Fig. 4. Analysis of pollen-tube growth in Ul and SI in semi-/n vivo pollinations of developmentally staged pistils. (A) The mean (±sd) 
style length (white bars) and the grand mean (±sd) of the 1 0 longest pollen tubes (black horizontal line) in Ul crosses of S. pennellii with 
cultivated tomato. (B) The mean (±sd) style length (white bars) and the grand mean (+sd) of the 10 longest pollen tubes (black horizontal 
line) in SI pollinations of S. pennellii with self pollen. Significant differences (P<0.05) within each pollination type (SI and Ul) for style length 
and for mean pollen-tube length based on Tukey's (honestly significant difference) post hoc tests are denoted within each figure by 
upper-case letters for style length and lower-case letters for mean pollen-tube length. The number of images analysed for each type of 
pollination is shown in parentheses. 



site of Ul pollen rejection. It should be noted that although 
expression of HT-A (in the upper stigma/style) is detected 
in S. lycopersicum (Fig. 7, arrowhead), this truncated HT-A 
lacks the asparagine-rich domain and thus is presumably 
non-functional in this fully compatible species (Kondo et ah, 
2W2a; Covey et al, 2010). 



Proteomic analysis of stigma/styles from key 
developmental stages 

iTRAQ profiling of -1 and -5 stigma/styles of SI acces- 
sion LA2560 resulted in the identification and quantification 
of a total of 2534 stylar proteins with a 95% CI from three 



Table 1. Results from four separate ANOVAs examining the effects of accession ID and bud stage on mean pollen-tube length and 
mean style length in SI and Ul crosses. 



Response variable 


Pollination type 


Factor 


df 


F value 


P value 


Mean pollen-tube length 


Ul 


Accession ID 


2 


0.077 


0.93 






Bud stage 


4 


8.78 


<0.000001* 




SI 


Accession ID 


2 


1.04 


0.31 






Bud stage 


4 


2.92 


0.03* 


Mean style length 


Ul 


Accession ID 


2 


1.72 


0.18 






Bud stage 


4 


72.84 


<0.000001* 




SI 


Accession ID 


2 


0.81 


0.37 






Bud stage 


4 


29.44 


<0.000001 * 



*Results of pairwise post hoc tests are shown in Fig. 4. 
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Fig. 5. SI proteins during pistil development in S. pennellii. (A) S. pennellii SI LA2560; (B) S. pennellii SC LA071 6. Immunoblots (1 .2 mg 
fresh weight/lane) show relative abundance and time of expression of S-RNase (S perv1 -RNase, an S-RNase previously identified in 
LA2560), HT-A and HT-B proteins in developing pistils. Protein sizes (kDa) are given on the left. Bottom panels show Coomassie- or 
silver-stained gels of the pistil extracts as loading controls. 




Fig. 6. Relative levels of S-RNases in S. pennellii accessions SI 
LA2560, SI LA1340 and SC LA0716. S-RNases in protein extracts 
(0.2 mg fresh weight/lane) from mature pistils were immunostained 
with anti-S-RNase C2 domain antibody. Silver staining of the 
blotted protein extracts acted as a loading control. 
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Fig. 7. Localized expression of HT-A protein in stigma/styles. 
Mature stigma/styles (lacking the ovary) of S. pennellii SI LA2560 
and SC LA071 6 and S. lycopersicum cultivar VF36 were divided into 
upper (UH) and lower (LH) halves, and proteins were extracted and 
separated on SDS/PAGE gels. Upper panel: immunostaining of pistil 
extracts using anti-HT-A antibody (1 .2mg fresh weight/lane). Protein 
size (kDa) is given on the left. Arrowhead on the right indicates the 
predicted truncated form of HT-A in cultivar VF36. Silver staining of 
the blotted protein extracts acted as a loading control. 
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Fig. 8. Biological process ontology classification of the top 1 00 differentially expressed proteins in inmmature (-5 stage) and in mature 
(-1 stage) stigma/styles. Proteins were categorized as having a signal peptide (SP) or not having a signal peptide (non-SP) based on 
N-terminal signal peptide prediction analysis using SignalP 4. Only three -5 SP proteins were identified and are not shown. Classification 
is according to the AmiGO gene ontology database (http://amigo.geneontology.org). 



separate iTRAQ experiments. From this group, 2015 proteins 
were identified with more than one peptide, allowing rela- 
tive quantification within each sample. For determination of 
the relative abundance of the proteins in each sample type, 
a selection parameter requiring the detection of any identi- 
fied protein in at least two of the three iTRAQ experiments 
was applied. A total of 1312 proteins met all three conditions 
(95% CI, more than one peptide detected, and detected in at 
least two of three experiments) and were used in quantita- 
tive analysis. Of these, 498 were found to be more abundant 
in stage -5, whereas 131 were more abundant in stage -1, 
based on a cutoff of 2-fold or greater difference. A summary 
of functional annotations of these proteins is shown in Fig. 8. 

The top 100 proteins that showed the greatest relative 
abundance at the -5 stage (immature stigma/styles) are 
shown in Supplementary Table SI. Prominent proteins in this 
group include histones (11%) and ribosomal proteins (58%), 
consistent with a phase of cell division and elevated cellular 
metabolism in immature stigma/styles. The top 100 proteins 
with a higher relative abundance at the -1 stage (mature 
stigma/styles) are listed in Supplementary Table S2. This 
group included proteins associated with lipid metabolism 
and transport, with 14% annotated as lipid transfer proteins. 
Accumulation of lipids is a hallmark of pistil maturation 
in most plant species with wet stigmas, including those in 
Solanaceae (Konar and Linskens, 1966; Dumas et al, 1978; 
Matsuzaki et al., 1983; Nakamura et al, 2003, 2009), Stigma/ 
styles at the -5 stage lacked a fully differentiated transmitting 



tract with its distinctive elongated cells, and also had very lit- 
tle or no lipid secreted onto the stigma surface, as determined 
from observations of sections of osmium tetroxide-stained 
stigma/styles (Fig. 9). In -1 and -2 stage stigma/styles, lipids 
are a prominent feature both on the stigma surface and 
between the cells of the now recognizable transmitting tissue 
(Fig. 9). These observations are consistent with the proteom- 
ics results and with the poor retention and growth of pollen 
on stigmas of stages younger than -3 in on-plant pollinations. 
Other proteins that were prominent at the -1 stage include 
arabinogalactan proteins, defensin, thaumatin, chitinase, 
HT-A proteins and S-RNases, among others (Supplementary 
Table S2). We performed a search for predicted N-terminal 
secretion signal peptides (SPs) using SignalP 4 with sequences 
of the most differentially expressed proteins at the -5 and -1 
developmental stages (SP; Supplementary Tables SI and S2, 
Fig. 8) using the most complete codons available for those 
sequences with a match in GenBank (www.ncbi.nlm.nih. 
gov/genbank/; 96/100 for the -5 stage and 95/100 for the -1 
stage). SP prediction analysis of differentially expressed pro- 
teins in the two stages revealed a major shift towards a secre- 
tory function in mature stigma/styles. Only 3% of the top 100 
proteins at the -5 stage stigma/styles were predicted to con- 
tain a SP, compared with 38% of proteins at the -1 stage. 

Known SI proteins were quantified in SI LA2560 stigma/ 
styles at the -1 and -5 stages. HT-A proteins were detected at 
higher levels at the -1 stage relative to the -5 stage, consistent 
with immunoblotting results (Fig. 5), and were represented 
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Fig. 9. Stigma/style development in S. pennellii. (A) Longitudinal sections of osmium tetroxide-stained, plastic-embedded pistils of S. 
pennellii SI LA2560 at -5 to -1 days before bud opening. An increase in stigma size and the amount of black, osmiophilic (presumably 
lipid) material in the stigma/styles is apparent during pistil development. The sections shown were taken from approximately the mid- 
sagittal planes of the stigmas. (B, C) Enlarged images of sections of the stigmas of -5 (B) and -1 (C) pistils, taken from the areas within 
the upper boxes in (A). There is little intercellular space among the stigma cells at -5, and only a small amount of osmiophilic material 
(arrowhead) is present in some cells at the surface. Papillae (P) are minimally developed. By -1 , the intercellular space among the cells 
of the stigma and upper style is greatly increased and filled with osmiophilic material which also covers the surface of papillae in the 
stigma (arrow). Globules of osmiophilic material are also abundant within papillae and in cells below the stigma surface (arrowheads). 
Papillae are well developed and have thickened cell walls. (D, E) Enlarged images of sections of stigma/styles at -5 (D) and -1 (E), 
taken from the regions in the lower boxes in (A). At -5 the cells in the central part of the style are relatively undifferentiated. By -1 the 
well-developed transmitting tract (TT) consists of markedly elongated cells. At both stages, osmiophilic globules interpreted as lipids are 
abundant within peripheral cells of the style. Scale bar in (B) also applies to (C-E). 



by at least four isoforms (Supplementary Table S3). HT-B 
protein was not detected in the proteomic analysis, suggest- 
ing a lower abundance of HT-B compared to HT-A. In addi- 
tion to the two S-RNases shown in Supplementary Table SI, 
peptides corresponding to six additional S-RNases were iden- 
tified in -1 stigma/styles of LA2560 (Supplementary Table 
S4), including one previously described in SI LA2560 (Pen-1; 



Covey et al, , 2010). In this group, one S-RNase was identified 
and quantified using a custom database generated from 454- 
based cDNA sequencing (RNA-Seq) data (TU1 56268; see 
Materials and Methods) and five were identified using a data- 
base collected from solanaceous S-RNase sequences anno- 
tated in GenBank. Since the latter contains only S-RNase 
sequences, search restriction parameters were increased to 
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compensate for potential bias, such that only those proteins 
identified with at least five different peptides were included in 
Supplementary Table S4. 

Discussion 

Developmental control of reproductive barriers 

Developmental regulation of reproductive barrier establish- 
ment is known to occur, since bud pollinations have been 
used to successfully overcome both interspecific and self- 
ing barriers in many species (Gradziel and Robinson, 1989, 
1991; Kuboyama et al, 1994; Sanchez et al, 2004; Singh 
et al, 2004). For example, UI barriers between cultivated 
tomato and either Solarium lycopersicoides or S. peruvianum 
were circumvented by bud pollinations using a stigma com- 
plementation procedure on otherwise non-receptive stigmas 
of pre-anthesis buds (Gradziel and Robinson, 1989, 1991). 
Bud pollinations have also been successful in overcoming 
reproductive barriers in other Solanaceae members, such 
as Nicotiana tabacum and Petunia hybrida (Sink and Power, 
1978; Kuboyama et al, 1994; Sanchez et al, 2004). 

Here we compare the timing of SI and UI barrier estab- 
lishment in pistils of three accessions of a single species, and 
present a detailed pollen-tube growth analysis during pistil 
development. The predictable sequential bud maturation on 
S. pennellii inflorescences and the development of a robust 
semi-/n vivo pollination system (Figs 1A, B and 2) contrib- 
uted to the success of our studies. The semi-z'n vivo system 
yielded consistent results by minimizing environmental and 
endogenous factors known to influence the results of on-plant 
bud pollinations in many plant systems, including wild spe- 
cies of tomato (Emerson, 1940; Linskens, 1975; William and 
Knox, 1982; Webb and Williams, 1988 and references therein; 
Gradziel and Robinson 1989). The variability of these factors 
can confound results and lead to conflicting inferences. For 
example, differences in temperature, availability of nutrients 
and photoperiod can contribute to the strength of SI reac- 
tions (Webb and Williams, 1988 and references therein). The 
semi-/« vivo pollination system allowed the reliable tracking 
of both SI and UI reproductive barriers during style develop- 
ment. Both types of barrier are absent 5 days before flower 
opening, and barriers are fully established by 2 days before 
flower opening (Figs 3 and 4). 

Our results as well as previous reports (e.g., Gradziel and 
Robinson, 1989, 1991; Kuboyama et al, 1994) show that the 
pistil becomes competent to support directional pollen-tube 
growth early during its development, but flower opening is 
delayed until the pistil establishes the molecular machinery 
that can discriminate between genetically compatible and 
incompatible pollen. Since S. lycopersicum pollen tubes grew 
within immature styles successfully to the ovary in S. pennel- 
lii pistils prior to barrier establishment, our results are more 
consistent with an active rejection (incompatibility) as the 
underlying mechanism for interspecific pollen rejection than 
a mismatch between pollen and pistil components (incongru- 
ency), with reproductive barriers layered on top of an other- 
wise compatible system. However, incongruency remains an 



explanation for failure of interspecific crosses in some sys- 
tems (Kermicle and Evans, 2005), and this phenomenon is 
likely to be widespread (Hogenboom, 1975). 

Sl-related proteins and the onset of UI 

Immunoblotting results (Fig. 5) show that the accumulation 
of HT-A protein coincides with the manifestation of UI bar- 
riers both in SI LA2560 and in SC LA0716. Proteomic analy- 
sis confirmed that HT-A protein was differentially expressed 
in the -1 stage compared to the -5 stage in SI LA2560 
(Supplementary Table SI). HT-B was not detected until later; 
a similar delay in the appearance of HT-B over HT-A tran- 
scripts in the stylar tissues of Solanum chacoense was reported 
by O'Brien et al, (2002). Although the HT-B antibody used 
is highly selective, it has a lower sensitivity than the HT-A 
antibody. Thus, the relative timing of HT-A versus HT-B 
expression is less clear than the overall correlation between 
HT-A expression and UI pollen rejection. Moreover, several 
lines of evidence point to a role for HT proteins in UI. First, 
the two tightly linked HT paralogs map to a pistil-side QTL 
controlling UI in the tomato clade (Bernacchi and Tanksley, 
1997; Covey et al., 2010). Second, HT expression patterns 
are correlated with UI at the species level. For example, all 
four of the SC red-fruited species (S. lycopersicum, Solanum 
pimpinellifolium, Solanum galapagense and Solanum chees- 
maniae) in the tomato clade lack UI barriers and lack both 
HT-A and HT-B proteins due to mutations in both genes that 
create early termination codons (Kondo et al, 2002a, 20026). 
The two SC green-fruited species (Solanum neorickii and 
Solanum chmielewskii) have UI barriers and express HT-A 
protein (Kondo et ah, 2002a). All tested SI accessions in the 
tomato clade have UI barriers and express both S-RNase and 
HT-A (Kondo et al, 2002a; Covey et al, 2010). The spatial 
correlation between UI and HT-A protein accumulation pre- 
sented here is also consistent with a role for HT proteins in 
UI (Fig. 7). 

Both S-RNase and asparagine-rich HT proteins are known 
to be required for SI because suppression of either S-RNase 
or HT protein prevents S-allele-specific self pollen rejection 
(Murfett et al, 1995; McClure et al, 1999; Hancock et al, 
2005; Puerta et al, 2009). In addition, some types of UI have 
been shown to require HT proteins in Nicotiana (Hancock 
et al, 2005). Our results showing significant accumulation 
of S-RNase prior to the onset of SI or UI in S. pennellii SI 
LA2560 (Fig. 5) are consistent with the regulation of the SI 
system in Nicotiana, where S-RNase transcript accumulation 
precedes HT transcript accumulation, and SI establishment 
requires sufficient levels of both proteins (McClure et al, 
1999). The difference between HT-protein accumulation in 
the upper versus lower half of the stigma/style (Fig. 7) corre- 
lates with the site of pollen rejection. The greater abundance 
of HT protein in the upper stigma/style is not likely to be 
attributable to higher expression in the stigma, since studies 
in Nicotiana showed moderately lower HT- transcript levels 
in stigma versus style (McClure et al, 1999). These results 
must not be interpreted to suggest that all active UI mecha- 
nisms are S-RNase-dependent. For example, SC S. pennellii 
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LA0716 lacks any detectable S-RNase protein or activity, but 
displays strong UI pollen rejection (Fig. 6, Liedl et al, 1996; 
Covey et al. , 2010). Thus, redundant S-RNase-dependent and 
S-RNase-independent mechanisms are likely to contribute to 
UI in the tomato clade, and the loss of SI (as is the case to SC 
S. pennellii LA0716) does not necessarily lead to a loss of UI. 

Changes in the stigma/stylar proteome associated with 
reproductive barrier establishment 

Proteomic analyses of developmentally staged SI LA2560 
stigma/styles revealed that the top 100 most abundant proteins 
at the earlier developmental stage (-5) featured histones and 
ribosomal proteins, reflecting a high rate of cell division and 
active metabolism (Supplementary Table SI, Fig. 8). At the 
later developmental stage (-1), when styles are able to reject 
incompatible pollen, the most abundant group of proteins 
is annotated as potentially associated with lipid metabolism 
[Supplementary Table S2, Fig. 8, e.g. lipid-transfer protein 
(LTP) and GDSL]. LTPs are proteins with a wide range 
of biological roles (Yeats and Rose, 2008) and these stylar 
LTPs may be functionally associated with the accumulation 
of large amounts of secreted lipid on the stigma and within 
the transmitting tract of stigma/styles immediately prior to 
flower opening (Fig. 9). Solanaceous species such as tomato, 
tobacco and petunia have wet stigmas with large amounts 
of lipid-rich exudates that peak during anthesis (Konar and 
Linskens, 1966; Dumas et al, 1978; Herrero and Dickinson, 
1980; Matsuzaki et al, 1983; Nakamura et al, 2003, 2009), 
and these exudates are critical for controlled pollen hydration 
(Wolters-Arts et al, 1998; Lush, 1999). 

Some of these LTPs, as well as other highly abundant pro- 
teins detected at this -1 stage, such as defensins and thau- 
matins, can be classified as cysteine -rich proteins (CRPs). 
Interestingly, some CRPs seem to have evolved functions in 
pollen-pistil interactions (Marshall et al, 2011). For exam- 
ple, Stigl is a stigma-specific CRP involved in exudate for- 
mation on stigmas and possibly in promoting pollen-tube 
growth (Tang et al, 2004; Verhoeven et al, 2005). Pistil LTPs 
facilitate pollen-tube adhesion and growth in the transmit- 
ting tract of lily and Arabidopsis styles (Park and Lord, 2003; 
Tung et al, 2005; Chae et al, 2007, 2009, 2010). CRPs that 
are more abundant in - 1 stage stigma/styles could play roles 
in defense (e.g. Lin et al, 2007; Hammami et al, 2009), since 
the stigma/styles at this stage are hydrated, nutrient-rich and 
exposed to the environment, and therefore accessible to path- 
ogens, or could function in pollen-pistil interactions (Kessler 
et al, 2010). Defensin-like proteins contribute to a number of 
reproductive processes, including pollen viability (Stotz et al, 
2009), pollen hydration (Muschietti et al, 1994), pollen- 
stigma recognition in sporophytic SI in Brassica (Kachroo 
et al, 2002), pollen-tube guidance to embryo sacs (Okuda 
et al, 2009) and the regulated release of sperm cells from pol- 
len tubes (Amien et al, 2010). At least 43% of proteins at the 
- 1 stage have predicted ER-targeting SPs, consistent with the 
importance of membrane-localized and secreted proteins in 
pollen-pistil interactions. Among these secreted proteins we 
noted a number of predicted cell-wall-remodelling proteins, 



including pectinases (four distinct putative pectate lyase 
proteins), and a polygalacturonase. Both of these classes of 
proteins presumably act to depolymerize the pectins in the 
primary cell wall and, perhaps more importantly, the middle 
lamella, a structure that may be less prominent in Solarium 
spp than in other taxa, although the intercellular material 
has been suggested to contain pectins (Cresti et al, 1976). In 
order for pollen tubes to penetrate the stylar tissues and reach 
the ovary, they need to pass between the transmitting tract 
cells. We suggest that in mature pistils these pectinases act to 
degrade the primary wall of the transmitting tract cells, or 
the pectin component of the intercellular material to facili- 
tate this process, as has been suggested previously for various 
pectinases (Marm-Rodriguez et al, 2002). A predicted expan- 
sin, another class of protein with a well-established role in 
cell-wall loosening, is also in the top differentially expressed 
100-1 stage proteins. Expansins in pollen have been previ- 
ously proposed to promote stylar penetration (Wu et al, 
2001), and their presence in styles could also promote stylar 
penetration by pollen tubes. 

Proteome characterization of SI LA2560 stigma/styles 
also confirmed the presence of S-RNases and HT proteins 
(Supplementary Tables SI, S2 and S3), which are known to 
contribute to SI in solanaceous species. Relative quantifica- 
tion confirmed that both of these classes of protein were more 
abundant in -1 stage stigma/styles than in -5 stigma/styles, 
correlating temporally with pistil maturity and the manifesta- 
tion of reproductive barriers. Our results therefore provide 
targets for future studies aimed at elucidating the multiple 
complex pathways that are Ul-associated. 



Supplementary material 

Supplementary material is available at JXB online. 

Supplementary Fig. SI. Style length during flower devel- 
opment. Fresh styles of -5, -4, -3, -2 and -1 stage flowers 
from S. pennellii accessions LA2560, LA1340 and LA0716 
were imaged immediately after placement of pistils on 
medium, and measurements were made using the ImageJ 1.33 
programme. N, number of images analysed. Error bars indi- 
cate standard error of the mean. 

Supplementary Table SI. The top 100 differentially 
expressed proteins in stigma/styles of S. pennellii SI LA2560 
at the -5 stage (immature stigma/styles). SP indicates 
N-terminal signal peptide prediction using the SignalP 4 pro- 
gramme; Y and N indicate secreted and non-secreted, respec- 
tively, according to the prediction. Ratio -1/-5 indicates the 
average relative quantification at the - 1 stage with respect to 
the -5 stage. SD is the standard deviation associated with the 
average ratio. Np indicates the number of peptides used for 
quantification, with a range indicated when different num- 
bers of peptides were used in different experiments. NR indi- 
cates number of iTRAQ experiments in which the protein 
was identified (out of three possible). 

Supplementary Table S2. The top 100 differentially 
expressed proteins in stigma/styles of S. pennellii SI LA2560 
at the - 1 stage (mature stigma/styles). SP indicates N-terminal 
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signal peptide prediction using the SignalP 4 programme; 
Y and N indicate secreted and non-secreted, respectively, 
according to the prediction. Ratio -1/-5 indicates the aver- 
age relative quantification at the -1 stage with respect to the 
-5 stage. SD is the standard deviation associated with the 
average ratio. Np indicates the number of peptides used for 
quantification, with a range indicated when different num- 
bers of peptides were used in different experiments. NR indi- 
cates number of iTRAQ experiments in which the protein 
was identified (out of three possible). 

Supplementary Table S3. HT-A peptides detected in stigma/ 
styles of S. pennellii SI LA2560. Peptide sequences represent 
HT-A proteins described in Covey et al. ,2010, as well as three 
additional HT-A peptides. Ratio -1/-5 indicates the average 
relative quantification in the -1 stage relative to the -5 stage. 
SD means standard deviation associated with the average 
calculation. NR indicates number of iTRAQ experiments 
in which the protein was identified (out of three possible). 
Np indicates the number of peptides used for quantification. 
NA means that the peptide was not used for quantification. 
Peptides are shown according to their position in the protein, 
from N- to C-terminus. The first four peptides are detected in 
S. lycopersicum cv. VF36, but not the remaining five peptides, 
which are distal to the premature termination site that pro- 
duces a truncated protein in this species. 

Supplementary Table S4. S-RNases identified in stigma/ 
styles of S. pennellii SI LA2560. The S-RNases in this table 
are in addition to the two S-RNases in Supplementary Table 
S2. -1/-5 indicates the average relative quantification in -1 
stage respect to the -5 stage. SD means standard deviation 
associated with the average calculation. Np indicates the num- 
ber of peptides used for quantification. NR indicates number 
of iTRAQ experiments in which the protein was identified 
(out of three possible). 
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